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Abstract: We numerically calculate the spatial correlators of the scalar and pseudoscalar 
operators F"^ and FF, in SU(3) Yang-Mills theory at zero and finite-temperature on the 
lattice. We compare the results over the distances < < ^ to the free-field prediction, 
to the operator-product expansion as well as to the strongly coupled large- A'^c = 4 super- 
Yang- Mills theory, where results are obtained by AdS/CFT methods. For Tc <T < 1.15Tc, 
both channels exhibit stronger spatial correlations than in the vacuum, and we give an 
explanation for this, using sum-rules and the operator-product expansion. The AdS/CFT 
calculation provides a semi-quantitatively successful description of the vacuum-subtracted 
F'^ correlator, renormalized in the 3-loop MS scheme, in the interval of temperatures 1.2 < 
T/Tc < 1.9, while the free-field prediction has the wrong sign. The FF and F"^ correlators 
are predicted to have the same functional form both at weak coupling and in the strongly 
coupled SYM theory. The Yang-Mills plasma does not meet that expectation below 2Tc. 
Instead we find that strong fluctuations of FF are present at least up to that temperature. 
We discuss the impact of our results on our understanding of the quark-gluon plasma. 
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1. Introduction 



Heavy ion collisions at RHIC have revealed properties of the quark gluon plasma that had 
not been widely anticipated (see [1] for an introduction). The ability of the produced matter 
to flow with little dissipation and to strongly quench energetic jets seemed to disfavor a 
description of the matter in terms of weakly interacting quarks and gluons. On the other 
hand, the constituent quark number scaling of the measured elliptic flow coefficient (see 
for instance [2] ) suggests that it is particles with the quantum numbers of quarks that are 
flowing in the expanding fireball. 

One of the central questions is thus whether the quark-gluon plasma at temperatures 
within reach of heavy-ion collisions is better described in a weak coupling expansion or 
whether a radically different computational scheme is more appropriate. The answer to 
the question could depend on the quantity, in which case it would be even more difficult 
to form a mental picture of the plasma. 

The strong elliptic flow and jet quenching observed in heavy ion collisions point to 
very strong interactions among the constituents of the plasma. Indeed the quantities most 
sensitive to interactions appear to be the dynamical ones, such as the shear viscosity ij 
in units of the entropy density s, which varies like aj^ between order unity and +oo 
with the coupling. Such dynamic properties of the plasma remain a challenge for lattice 
calculations (see [3] for a review). In the shear channel the most accessible transport 
property is dijop{u)/u with of order T, where p is the spectral density. For a weakly 
coupled system obeying the /-sum rule [4], this provides a measure of the mean square 
velocity of the quasiparticles responsible for the transverse transport of momentum. A 
value much below unity would rule out the possibility of these quasiparticles being light 
quarks or gluons. 

Static quantities on the other hand, while often providing a less clear-cut test of the 
importance of interactions, are directly accessible in the Euclidean formulation of the the- 
ory. The thermodynamic potentials, for instance, have remained a challenge for pertur- 

bative methods [5, 6], even though certain resummation schemes lead to more stable pre- 
dictions [7, 8]. The example of strongly coupled = 4 super- Yang- Mills (SYM) theory, 
where the entropy density is only reduced by a factor 3/4 with respect to the non- interacting 
case [9] , shows that only a highly accurate agreement of the weak coupling expansion and 
non-perturbative lattice data can warrant the conclusion that the plasma is dominated by 
weakly coupled quark and gluon quasiparticles. 

Other static quantities on the other hand appear to be quite well described by weak 
coupling techniques. A convincing example is the spatial string tension, for which the di- 
mensional reduction program works well [10, 11]. As another example, the fluctuations of 
quark numbers [12] appear to approach remarkably early the Stefan-Boltzmann limit. Re- 
cently the expectation values of other twist-two operators (besides the energy-momentum 
tensor) have been proposed [13] as diagnostic tools for the effective strength of interactions. 

In this paper we calculate non-perturbatively spatial correlators of two dimension-four 
operators, the trace anomaly 6{x) and the topological charge density q{x) in the SU(3) 
gauge theory. The range of distances covered by the calculation is < r < ^ . These are 
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also static quantities that are directly calculable on the lattice. Furthermore, quite a lot is 
known about these correlators in the QCD vacuum, going back to the original QCD sum 
rules studies [14, 15]. And thirdly, they are computable in the large- A'ci strongly coupled 
A/" = 4 SYM theory by AdS/CFT methods [16, 17, 18]. Thus we have the possibility to 
compare the lattice data to two 'caricatures' of the plasma, one being non-interacting gluons 
and the other being a very strongly coupled non-Abelian plasma. As we shall see, once 
the vacuum contribution has been subtracted these two caricatures lead to qualitatively 
different predictions for the correlators. 

Parallel to the question of the weak- or strong-coupling nature of the quark-gluon 
plasma, lies the question of how similar non-Abelian relativistic plasmas are. This consti- 
tutes a very interesting question in itself. In addition, the possiblity to compute real-time 
quantities in strongly coupled theories amenable to AdS/CFT computations and "port" 
them to QCD provides a strong phenomenological motivation. The best known examples 
of this strategy are the shear viscosity to entropy density ratio [19] and the jet quenching 
parameter q calculations [20]. As evidence in favor of the strategy, the authors of [21] 
conclude that the overall agreement of the screening spectra of QCD and the J\f = 4 SYM 
theory is rather good, although the low-lying screening masses are overall a factor 1.9 or so 
larger in the strongly coupled SYM theory. They therefore suggest that the QCD plasma 
around 2Tc is most similar to the = 4 SYM plasma at an intermediate value of the 
't Hooft coupling A. In order to find the effective coupling which leads to the best match 
(defined by a set of physical quantities) between the theories therefore requires knowing the 
properties of the A/" = 4 SYM plasma at intermediate values of the coupling, presumably 
as hard a problem as determining those of the QCD plasma. However, in the SYM theory 
one has the advantage of being able to expand the observables in A and in 1/A, opening 
the possibility to interpolate to intermediate couplings [21]. 

The outline of this paper is as follows. In section (2) we define the relevant operators 
and their correlators, and give the basic free-field theoretic predictions. Section (3) contains 
the AdS/CFT calculation of the same correlators in the strongly coupled J\f = 4 SYM 
theory. In section (4) we describe the lattice calculation of these correlators, including 
a new way to normalize the topological charge density for on-shell correlation functions. 
The results are compared to weak- and strong-coupling theoretical predictions. Section (5) 
discusses what values of the 't Hooft coupling best match the gluonic and the SYM plasma. 
We finish with a summary of the lessons learnt and an outlook in section (6). 

2. Definitions and theoretical predictions 

In this section and in the following, we use Euclidean conventions, since the calculation of 
correlators will be performed on the lattice. We consider the SU(iVc) gauge theory without 
matter fields. 



We focus on two operators in this paper. The first is the (anomalous) trace of the energy- 




(2.1) 
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momentum tensor, 

0(x)^r^^(x) = ^F;,F;,, ^(5) = -6/ + ..., 6=3^. (2.2) 

The second operator is the topological charge density. It is defined as 

= ^e^upaHF^Ax)F,,{x)} = (2.3) 

where F^i, = gF^j^t"-, ii{t°'t'^} = —\5ab and F^jy{x) = ^efj,,^prjFp^{x). The normalization is 
chosen such that the value oi Q = J d'^x q{x) on a self-dual configuration is an integer. For 
later use we also introduce the operator 

9oo{x) = ^F-^F-^-^FoiFoi. (2.4) 

In the thermodynamic limit, (^oo)t-o = e + p while {9)t-o = e — 3p, where e and p are the 
energy density and pressure respectively and the subscript T — means that the difference 
of the thermal expectation value and the vacuum expectation value is taken. For O = d or 
q, we will consider the static connected correlators at finite temperature T = 1/Lo, 

Coo{r,T) ^ (O(0,r)O(0))c = ^(^^^ k''"'' ~ [e-^°^a(0)]) 

(2.5) 

Often, to emphasize the thermal effects on the correlator, we will subtract the zero- 
temperature correlator. 



Gooir, T) = Coo{r, T) - Coo{r, 0) . (2.6) 

If one expresses the traces of Eq. (2.5) in a basis of energy eigenstates, this subtraction has 
the effect of removing the vacuum contribution. 



2.1 Short-distance behavior 

In this section we review the available perturbative results for the correlators (2.5) as well 
as our knowledge of their long-distance behavior. 



2.1.1 Zero temperature 

The two-point functions of the trace anomaly and the topological charge density are to 
leading order 

(87r6a,)-2(0(x)^(O))noop = - (-)' (^(x) g(0))noop = , (2-7) 

\as / TT [X ) 

where cIa = — 1 is the number of gluons. The correlators were calculated to two-loop 
order in [22], but we will not exploit that result. 
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2.1.2 Finite temperature 

The two-point functions of the trace anomaly and the topological charge density to leading 
order read [23] 

/( \ tXj T Iwf 1 



m,neZ V H N 



(2.8) 

where = x + nLQCQ. Their operator-product expansion (OPE) reads [14, 23] 



i8.ba.r^{eix)em = - (g)^,(x) ,(0)) + -- - = ^- ^^-i^ + -- -(2.9) 

and the dots refer to O(^) terms. In fact, the OPE of the two correlation functions 
appearing in Eq. (2.9), treating the Wilson coefficients to leading order in perturbation 
theory, differ only at 0{r^) if one restricts the terms on the right-hand side to operators 
whose vacuum expectation value does not vanish [15]. Since the term cancels exactly 
when taking the difference of two temperatures and {0oo)t-o > 0, (^)t-o > Oj Eq. (2.9) 
implies that the gluon plasma is always more screening than the vacuum of the theory at 
sufficiently small distances. 

2.1.3 Spectral functions 

The free spectral functions for the trace anomaly and the topological charge density are [24] 

(87r6«,)-V,(a;,g,r) = - (g)%^^(a;, g, T) = ^ (a;^ - g')' X([l], a;, g, T) , 
-^.r n N 2r , sinh(a;-Fg)/4r 

XI, a;, q, T) = e{q - u) + — log — . 2.10 

q q smh|a; — g|/4r 

The spectral functions are related to the static correlators by Fourier transformation, 

(O(0,r)C»(0)) = lim f -^e''!" H dwe-'^ ^^"^'"^'li • (2.11) 
\ ^ ' ^ ^ ^/ e^oj (27r)3 Jq tanha;/2r ^ ' 

The parameter e serves to regulate the integral over frequencies at large to, which would 
otherwise diverge. We will come back to this regulator in section 4. 

2.2 Long-distance behavior 

At long distances, the vacuum correlators of 6 and q are dominated by the scalar and 
pseudoscalar glueballs respectively. The most recent lattice results for their masses are 
roMo++ = 3.96(5) [25] or 4.16(11)(4) [26] and roMo-+ = 5.93(16) [27] or 6.25(6)(6) [26], 
where ro ^ 0.5fm is the Sommer reference scale [28]. The coupling of these states to the 
local operators 6 and q, s = (vac|0|O++) and p = {Yac\q\0 ^), have also been calculated 
recently [26, 25]. 

The screening masses, which determine the asymptotic exponential fall-off of the finite- 
T correlators, are also known to some extent. The operators 9 and q belong to irreducible 
representations (irreps) of the S0(3) rotation group, x parity and x charge conjugation. 
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At finite temperature, the symmetry group of a 'z-sUze' (for states at rest, Px = Py = 
and given Wn, the discrete momentum in the direction of length l/T) is reduced to 
RxSO(2)xP2 X C, where R is the EucHdean-time reflection and P2 is the reflection inside 
an xy plane {{x,y) — ^ {x,—y)). In general, an operator forming an irrep of the zero- 
temperature theory gets decomposed into several irreps of this reduced symmetry group. 
In our case, 9 and q simply become the scalar and pseudoscalar representations of the z-slice 
symmetry group. The former is invariant under all the symmetries of the z-slicc; the latter 
is too, except for being odd under the R and P2 operations. On the lattice, these irreps are 
generically further reduced to crystallographic irreps. In our case they are labelled A'^'^ and 
Ai'^ [29]. A recent result for the mass gap in the scalar sector is m^++{T)/T = 2.62(16), 
2.83(16) and 2.88(10) respectively at 1.24rc, LeST^ and 2.20re [23]. Datta and Gupta find 
m^-+{T)/T = 6.32(15) both at about 1.5Tc and 2.0Tc [29]. So the asymptotic screening is 
much stronger in the pseudoscalar sector than in the scalar sector. This ordering persists 
at all temperatures according to a recent next-to-leading order perturbative analysis, in 
spite of a change in the nature of the lightest scalar state [30] . 

3. AdS/CFT Results for = 4 Plasma 

We now turn to the calculation of spatial correlators in a maximally supersymmetric 
strongly coupled plasma using gauge-gravity duality [16, 17, 18]. We will find that the 
thermal correlators in the J\f = 4 plasma are identical for the operators and FF, as 
each of these operators is dual to a simple massless scalar field. It is interesting that the two 
correlators also coincide at weak coupling, where they are given by a two-gluon exchange 
diagram, and therefore coincide with the pure Yang-Mills result (2.8). 

These correlators have been previously studied in momentum space in [31]. See also [32, 
4] for discussion of finite temperature stress tensor and R-current correlators in momentum 
space. An outline of the computation is: 

1. letting O denote either F^ or FF, we note that the field theory operator O is dual 
to a massless bulk scalar field cf). For F^ this field 4> is exactly the type IIB dilaton, 
whereas for FF (f) it is the Ramond-Ramond axion Cq. 

2. We then compute the spectral density piuj, k) for the operator O using finite-temperature 
AdS/CFT. This involves numerically solving the bulk equations of motion for ^ in a 
black brane geometry. 

3. Finally, we Fourier transform this spectral density to obtain the Euclidean correlator 
in position space. 

Each of these steps is explained in more detail below. Throughout we will expand fields 
on each constant-radius slice in Fourier space, ^ ~ ^(r)e~*'^*+*'^^. We take the spatial 
momentum to be in the z direction. With an eye on numerical evaluation, we will often 
work with dimensionless momenta and positions, which we denote with an overbar: 



The relevant black brane metric for J\f 
written 



4 SYM at finite temperature on M^'^ can be 



2„2 



ds^ = {nRTYr 
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dt^ + dxl 



dr'^B? 



+ 



1 - 



(3.2) 



where R is the radius of the bulk AdS space and T the temperature of the black brane, 
with the horizon at r = 1 and the AdS boundary at r — > oo. 

3.1 Flow Equation and Numerical Evaluation 

We now let O denote either or FF. In both cases the relevant bulk action for the field 
dual to O is simply that of a massless scalar 



S= 

2a 



]-j d^x^g{Vct>f, 



(3.3) 



For these operators supersymmctry guarantees that the vacuum two-point function is inde- 
pendent of the coupling [31], and thus the normalization a can be conveniently determined 
by demanding that this correlator as computed from gravity agrees with the free-field 
expressions (2.7). We find for both 

1 1 Ar2 



2a r 



(3.4) 



2aF2 ^^pp 

The spectral density p is proportional to the imaginary part of the finite temperature 
retarded correlator, p = — ^Im(G_R). An extensive literature exists on the evaluation of 
this quantity from AdS/CFT [33, 34, 35, 36]. Wc will use the flow formalism developed in 
[37], which we briefly review here: consider the function x(r, k), defined as 



II /■. k) 



U{r,k) = —,/^g''drct>. 



(3.5) 



Here Il{r,k) is the momentum conjugate to the bulk field 4>{r,k). The bulk equation of 
motion for implies that the x(r, k) obeys (on any metric) the first-order flow equation 

,2 / ^2gzz^ 



rX = 




X y7 



1 - 



Ijj2gtt 



Sri, = — 

a 



-9 

drrgtt 



(3.6) 



Furthermore it follows from real-time AdS/CFT [38] that the retarded correlator Gr in 
the dual field theory is obtained from the boundary value of x(r, k): 



Gnik) = — lim iLox{r, k) 



Im[GR{k)] = - lim ujRe[x{r,k)]. (3.7) 



The initial conditions at the horizon r = 1 are fixed by the infalling wave condition to be 
X(r = 1) = S^(r = 1). 

We now plug in the metric (3.2) and define a dimensionless function xhy x= ^J^^s ■ 
We obtain the flow equation 

drX 



2iuj 



2(1 



x' 



W ( 1 
1-3^ (1- — 



(3.8) 



This equation must now be integrated from x{z = 1) = 1 to the AdS boundary at z = oo, 
where it determines the AdS/CFT response. Some technical details on the numerical 
integration are given in Appendix A. 
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3.2 Fourier Transform 



To obtain a Euclidean correlator from the spectral density, we use the identity [39]. 

roc p izt coshfcjfr— f 

Gs(0; r, x) = / —- p{u, k) (3.9) 

Jo J {^T^r sinh(^) 

Assuming rotational symmetry in the spatial directions (i.e. p{u, k) = p{uj, \k\)) to perform 
the angular integral and switching to dimensionless momenta, we obtain 

GE(0;f,|x|) = 87r2rM dw dlfel p(u;, fc) sin(|fc||a;|)- ' " (3.10) 

Jo Jo X smh(a;7r) 

We note at this point that we are primarily interested in equal time correlators, i.e. those 
for which r is in the equation above. However, at small time separations the nonzero value 
of r acts like a UV cutoff on high-frequency modes, suppressing them as e~'^'^. To achieve 
arbitrarily small r we would need to know p at arbitrarily high cj, whereas numerically we 
are necessarily limited to finite uj^. Since the position-space, equal-time correlator is finite, 
one could repeat the calculation with several r values and extrapolate to r = 0. Here we 
however restrict ourselves to using a small, finite 'regulator' r ^ r. Since our goal is to 
compare the correlators to those computed on the lattice, where the region of very small r 
is in any case problematic due to discretization errors, this approach will prove sufficient. 

As a check on the numerical Fourier transforms themselves, we compute the correspon- 
ding Fourier transform in the free theory at finite temperature starting from the analytic 
expression for the spectral density Eq. (2.10); in this case an analytic expression also exists 
directly in position space, Eq. (2.8), providing a non-trivial check on the accuracy of the 
Fourier transform. In both cases we subtract the zero-temperature contribution, which, as 
mentioned above, is independent of the coupling. The step sizes in O and k are both 0.1, 
and the chosen range of integration is to 20. The result of this numerical integration is 
shown in Fig. 8, where we have fixed r = With x ranging over values much larger 

than r, the Figure shows the departure of the Fourier-transformed correlator from the di- 
rect evaluation of the x-space expression (2.8). The achieved accuracy is sufficient for our 
purposes. The largest discrepancy occurs at short distances, where the sensitivity to the 
discretization step and the finite range of uj is greatest. 

To illustrate the dependence on the regulator r, we compare the correlator Ce)e)(r, r, T) — 
Cc)c)(r, r, 0) for r = l/27rT with Gooi^^T) in the free case (Fig. 9). We see that for 
r > 5/27rr, the correlators coincide to the accuracy needed for our purposes. 



4. Correlators on the lattice 



In this section we describe the lattice setup and the numerical results obtained by Monte- 

^Note this is an advantage to using the real-time formahsm described here, as there would be no such 

exponential suppression of high frequency modes if we were to compute the position space correlator by 
summing the Euclidean momentum space correlator over Matsubara modes. 



-7- 



Carlo simulations. We employ the Wilson action [40], 

Ss = \j2 tr{l-P^.(x)}, (4.1) 

where the 'plaquctte' P^i^ is the product of four link variables Un{x) around an elementary 
cell in the {fj,, v) plane. 

As a simulation algorithm, we use the standard combination of heatbath and over- 
relaxation [41, 42, 43, 44] sweeps for the update in a ratio increasing from 3 to 5 as the 
lattice spacing is decreased. The overall number of sweeps between measurements was 
typically between 4 to 12. 

4.1 Discretization and normalization of the operators 

A choice has to be made for the discretization of the operators 6{x) and q{x). Here we use 
the specific discretization 

OLix) = -Xs{9o)^^^ ^^Retr[i?,,(x)F^,(x)] , (4.2) 

Ql{x) = 'If^f' e/xz^pcxtr {x)Fp^ {x)] , (4.3) 

where the (antihermitian) 'clover' discretization of the field-strength tensor Ffj_^(x) is de- 
fined in terms of the link variables in [45]. In this work we feed in 'HYP smeared' link 
variables [46] into the definition of F^i^^x). The name stems from the fact that the elemen- 
tary link variable is replaced by an average of Wilson lines which remain in the adjacent 
elementary hypercubes. We kept the original parameters [46] and used the projection onto 
the SU(3) group of the Wilson- line average described in [47] . 

At the quantum level, the normalization of these lattice operators differs from the naive 
normalization. Indeed, even though the anomalous dimension of these operators vanishes, 
a finite renormalization of the operators survives, which has to be compensated for in order 
to ensure that their on-shell matrix elements approach their continuum limit with O(a^) 
corrections. 

J -2 

In Eq. (4.2), is the lattice beta-function which describes by how much the lattice 

spacing shrinks when the bare coupling is reduced. While asymptotically it is governed 
by the first two universal beta-function coefficients, we work in the region ^Iq ~ 1 and 
therefore employ a non-perturbatively determined beta-function. Specifically we use the 
quantity ro/a (the Sommer reference scale) as a function of gg , as computed in [48] and 
parametrized in the appendix of [49]. By taking one derivative of the parametrization, we 
obtain ^j^. 

a log a 

The trace anomaly in our chosen discretization still requires the additional normaliza- 
tion factor Xsido)- The latter is fixed by calibrating against the 'canonical' discretization of 
9{x). This discretization arises from differentiating the lattice action (4.1) with respect 
to the bare coupling. 



2 



d\og( 
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By requiring that e — 3p he independent of tlie discretization, i.e. {9l{x))t — {0Lix))o be 
equal to {9'j^{x))t — we determine Xs(5'o)- Here we choose A^,- = L^/a = 6 to do 

this matching. The results are well parametrized by 

3257 

>"fa°>= l-0.7659,§ - 5-9»£''£6.72. (4.5) 

The error varies from 0.004 at /? = 6.0 to 0.008 at /? = 6.72. We have not investigated 
systematically the dependence of Xs on the value of Nt used for the matching. This un- 
certainty is not included in the above error estimates, and requires repeating the matching 
calculation at A^^- = 8 or 12. For our purposes in this work, this uncertainty will not prevent 
us from drawing conclusions when comparing the trace anomaly correlator to theoretical 
predictions. 

4.1.1 Normalization of the topological charge density 

The procedure we use to normalize qiix) is somewhat new to our knowledge, and therefore 
we describe it in some detail. We again exploit the fact that there exists a discretization 
q']^{x) for which the normalization is known exactly. This is the definition based on the 
overlap operator [50] which satisfies the Ginsparg- Wilson relation [51]. Indeed, Q'l{x) = 
^x^ii^) ^^^T^ guaranteed to be an integer on every gauge field configuration, because 
it counts the difference of the number of right-handed and left-handed zero-modes of the 
Dirac operator [52]. 

We normalize our discretization of q{x) by matching the value of its vacuum two-point 
function with the same correlator obtained with the overlap-fermion-based discretization 
of q{x). The latter correlator was obtained in [53], and we use the numerical data of that 
article to fix the normalization of our discretization. More precisely, the quantity matched 
is Cqq (r, 0) ; this removes the largest part of the uncertainty in the lattice spacing in 
physical units. 

Specifically, we choose the matching distance to be f/ro 0.68, and use the data of [53] 
on the finest lattice (£^3 ensemble at a = 0.082fm). The matching distance f was chosen 
such that on this lattice, f/a = 3\/2. The lattice spacing in [53] is specified by the string 
tension, and wc have used the factor toa/o- = 1.1611(95) (based on the compilation [25]) 
to convert physical distances in units of tq. 

Given our normalization strategy, it is convenient to split the normalization factor of 
the operator into two parts: 

ZgiP) = Z,(/3,ef) • xiP)- (4.6) 

where we choose Pref = 6.2822. The advantage of this separation is that only Zq{l3rei) 
depends on the overlap-based correlator. The latter has been calculated to much lower 
statistics than our computationally cheaper correlator. We find 

^,(/3,ef) = 1.55(8), (4.7) 

where the uncertainty comes almost entirely from the overlap data. We then obtain, by 
matching the correlator at other values of (3, 

X(5.903) = 1.25(4) Wtch/ro = 0.94 (4.8) 
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X(6.018) = 1.125(24) Wtch/ro = 0.77 (4.9) 

X(6.200) = 1.030(15) wtch/r-o = 0.68 (4.10) 

X(6.408) = 0.959(20). r„,atch/ro = 0.60 (4.11) 

X(6.720) = 0.958(50). r^atch/ro = 0.60 . (4.12) 

When necessary, we use linear interpolations in r of the function r^Cqq{r,0) to match 
different lattice spacings. These % factors are well parametrized as 

0.8112 - 0.73885^ o r, 2 no^ 

= P = 6/on • (4.13) 

^^^^ 1-0.9364C/2 ' ^ /yo \ J 

The absolute error is about 0.041 at /3 = 5.9, 0.014 at (3 = 6.1, 0.010 at /? = 6.4 and 0.026 
at /? = 6.7. 

4.2 The vacuum correlators (T = 0) 

The vacuum correlators of the trace anomaly and the topological charge density (multiplied 
by r^) are displayed in Fig. 1. The overall normalization uncertainty coming from Xsido) 
and Zq{(3) are not included in the error bars on the picture; they are given above. Only 
data points with r/a > 4 are displayed. We have measured the correlator along the lattice 
axes (1,0,0), (1,1,0) and (1,1,1) and averaged the results over all equivalent directions. 
Some of the raw data is given in Tables (1) and (2). 

The plots also show the perturbative prediction (2.7) at small distances. The three 
lines give an indication of the renormalization scale uncertainty: they correspond to ^ = 
|:, ^, ^. For that purpose we have used the resuh A^=°Vo = 0.602(48) [54]. Based on 
these figures, it is very plausible that our data agrees with perturbation theory at short 
distances, but data at smaller lattice spacing is needed for a more stringent test. 

Such vacuum correlators are of course interesting in their own right. In particular, they 
can be used to test models for low-energy QCD, such as the instanton- liquid model [55] 
or the more recent holographic models of hadrons [56]. See [57, 58] for model calculations 
of gluonic vacuum correlators. A detailed comparison of the latter with lattice data will 
be carried out elsewhere; for the time being we simply note that around r = 0.4fm the 
lattice data points lie on a convex curve in the scalar case, and a concave curve in the 
pseudoscalar case. This observation is qualitatively consistent with the instanton-model 
calculations of [57], where it is interpreted as an attraction/ repulsion respectively in the 
scalar and pseudoscalar channels. 

4.3 Finite-temperature correlators 

The finite-temperature correlators of the trace anomaly minus their zero-temperature coun- 
terpart are displayed in Fig. (2). Partial results for this quantity were already published 
in [23] . A sample of the raw data is presented in Tables (3) and (4) . 

Starting from the higher temperatures (bottom panel), we note that the lattice cor- 
relators arc gradually approaching the free-theory prediction (Eq. 2.7) as the tempera- 
ture is raised, as one expects on the basis of asymptotic freedom. In order to display 
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Eq. (2.7) on the figure, we have set the renormahzation scale to = 7r(T + i) and used 
a1^~'^Vo = 0.602(48) [54]. However, at the temperatures where simulation data arc avail- 
able, the subtracted correlator is negative at all separations r - unlike the free correlator. 
This means that the plasma screens the fluctuations of the operator 9 more than the 
vacuum does. 

We have checked for discretization errors by calculating Gg0{r, 1.24Tc) at several lattice 
spacings. This means that N-j- = {aT)~^ is varied and (3 tuned so that the temperature is 
kept fixed. The result is shown in Fig. (5), top panel. To a good approximation, the data 
fall on a single curve. It is not presently our intention to carry out a systematic continuum 
extrapolation of G0e{r,T) in l/N^- Rather we want to show convincing evidence that 
the conclusions we draw from finite lattice spacing data are not affected by discretization 
errors. 

In the figures, we have not included the statistical uncertainty on the non-perturbative 
normalization factors. However the latter cancel in the ratio 

which is displayed in Fig. 3, top panel. We have multiplied this quantity by cIa so as to 
give it a finite limit when Nc ^ co and by (Tr)"^, so that the expected short-distance 
behavior is ~ "^(m)) where = 0(l/r). The graph shows that GQ0{r,T) falls off like 
as expected from the OPE, for 3 < 2TrTr < 5. Beyond this interval, it falls off faster to 
zero. 

The on-shell correlation functions of 6 are renormalization-group invariant. However, 
if we want to compare the result with the = 4 SYM correlator of F'^ described in 
section (3), we have to divide out the factor of P{g)/2g that multiplies the Yang-Mills 
operator 9. Since the beta-function is scheme-dependent, this means that the renormalized 

correlator in the Yang-Mills theory is scheme-dependent, too. We choose the three- 
loop MS scheme, and evaluate the coupling ag at the scale fi{r,T) = |7r(r-|- 1/r). A few 
remarks may be useful to motivate this choice. At r = 1/T, this corresponds to /j, = 
it makes the one-loop prediction for the 9 correlator approximately go through the lattice 
data at T = 0, sec Fig. (1). At short distances, we expect r to provide the harder scale and 
therefore the appropriate // to be dominated by r. For r > T, we expect the momenta of 
the two exchanged gluons to be of order ttT. The chosen expression for n is then a simple 
interpolation between these two regimes. 

We thus obtain the lattice data for the F^ correlator in the MS scheme, see the bottom 
panel of Fig. 3. We are then in a position to carry out a parameter-free comparison with 
the one-loop result Eq. (2.7) and the AdS/CFT result. In the range of temperatures 
1.2 < T/Tc < 1.9, the lattice correlators are in semi-quantitative agreement with the 
corresponding correlator calculated in the strongly coupled A/" = 4 SYM theory. The 
lattice data are negative at all 1/2T < r < 1/T, and this is in contrast with the free-field 
theory result, which is positive in that range. The data at T > 2Tc however does suggest 
that the non-perturbative correlator gradually moves towards the one-loop result as the 
temperature increases, as expected from asymptotic freedom. 
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Coming back to Fig. 2 (top panel), the thermal fluctuations become stronger as the 
temperature is lowered toward Tg, to the point where the subtracted correlator is positive 
over a wide range of distances r. Unsurprisingly this fact is accounted for neither by the 
weak coupling predictions, nor by the conformal M = 4 SYM result. The correlation 
between the fluctuations of 9 is strongest at l.OlTc, and drops again as one moves away 
from the transition below Tc- The interpretation of the data is helped by the fact that the 
overall fluctuations of 9 are related to thermodynamic properties via the sum rule [59, 60] 

J d'x {9{x)9{0)Unn,T-O = (4-15) 

Since we know that (e — 3p)/T^ rises very steeply between Tc and l.lTc [61], Eq. (4.15) 
indicates that the fluctuations of / d'^x 9 are strongest in that range of temperatures. Since 
the Wilson coefficients in the OPE are negative, there is a negative contribution to the 
LHS of Eq. (4.15) from the short-distance part of the correlator. Therefore there has to be 
an enhancement in the 9 two-point function at intermediate or long distances in order to 
account for the positive sign of the RHS (note however that here we restrict ourselves to 
equal-time correlators). On the other hand, above l.lSTc, where the RHS of Eq. (4.15) is 
negative, there is no necessity for Gqq to be positive at any non-vanishing separation. Our 
data shows that indeed Goo at intermediate distances r ~ 1/T is negative for all available 
temperatures above 1.2Tc. This differs from the correlator Gee of the energy density [23], 
which is positive at intermediate separations. 

4.3.1 Topological charge density correlator 

We now discuss the topological charge density correlator, starting from the high temper- 
ature end. An example of Gqq at high-temperatures is displayed in Fig. (6). We see that 
qualitatively, the correlator, in the interval where data is available, resembles its free-theory 
counterpart: —Gqq is negative at short-distance, as predicted by the OPE, and then pos- 
itive at intermediate distances. On the basis of the pscudoscalar screening masses, we 
expect Gqq to approach zero from below, and there is a hint at 3.30Tc that indeed it does. 

If we now lower the temperature, as shown in Fig. (4, bottom panel), we see that the 
range of distances where —Gqq is positive grows, and that its maximum value also grows. 
The OPE tells us that at sufficiently short distances, —Gqq must be negative. Below 
1.6rc its maximum is no longer visible in the data; we conjecture that it is located at too 
short distance r for us to see it in the lattice data (at short distances, we are limited by 
discretization errors). 

As we lower the temperature further (top panel of Fig. (4)), —Gqq{r,T) at fixed r 
continues to grow. It hits a maximum between 1.02 and l.OGTc. Thus similarly to the 
trace anomaly, the topological charge density exhibits strong spatial correlations near Tc. 
How strong they are is better illustrated by taking the ratio of the finitc-tcmpcrature to 
the zero-temperature data, Fig. (7). Here one clearly sees that for r of order 1/T, the 
spatial correlation of topological charge density fiuctuations is about twice as strong near 
Tc as in the vacuum. A technical advantage of this ratio is that the overall normalization 
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cancels out, and secondly that we expect a partial cancellation of the discretization errors 
to take place. 

To check that these large correlations are not a cutoff effect - there is after all a large 
cancellation taking place at short distance in the subtracted correlator - we have repeated 
the 1.24Tc calculation at a finer lattice spacing. The comparison is shown in Fig. (5). 
We see that the different data sets fall on top of eachother within errors in the interval 
0.4 < Tr < 1.2. We thus conclude that the strong, finite-temperature induced enhancement 
of the correlation is a true physical effect. 

There are significant differences between the scalar and the pseudoscalar channels in 
the lattice data at distances 1/2T < r < T. This is unlike the strongly coupled M = 4 
SYM theory, where the and FF correlators are identical. It is also unlike the free field 
prediction. In the OPE framework, this difference requires operators of dimension 6 or 
higher to overwhelm the expectation value of the stress-energy tensor in Eq. (2.9)^. This 
would in turn imply the breakdown of the OPE as an asymptotic expansion. 

Our discovery of large spatial correlations in the topological charge density fluctuations 
in the vicinity of Tc is qualitatively in agreement with the results of [62] . The latter showed 
a strong suppression of the topological susceptibility just above Tc, using a method based 
on the semi-classical identification of topological charges. This suppression is particularly 
dramatic at larger Nc values, but even for SU(3) it amounts to a factor of about 0.54(4). 
This implies that — / d'^x{q{x)q{0))T-o > 0, and therefore there has to be a range of 
separations x where —{q{x)q{0))T-o is positive; this is what we are seeing in the data. Note 
that the short-distance singularity of {q{x)q{0))T-o oc a'^/x^ gives a finite contribution 
when integrated over space-time. This is in contrast with the topological susceptibility 
itself, Xt = f d'^x{Q{x)q{0))o, which has to be defined with care [63] if it is to remain finite 
when the cutoff is removed. 

5. The effective coupling in the plasma 

We have found that the strongly coupled SYM theory has an F'^ correlator similar to the 
pure Yang-Mills theory in the deconfined phase below 2Tc. To summarize the procedure, we 
have calculated the 9 correlator on the lattice, which contains a factor {(3{g)/2g)'^ relative 
to the F'^ correlator. This factor makes it renormalization-group invariant in the pure 
Yang-Mills theory. We used the 3-loop MS scheme for the beta-function to convert the 
lattice 9 correlator to the F^ correlator, and found semi-quantitative agreement between 
the theories in a range of temperatures. For r = 1/T, the values for our chosen running 
coupling are 

as{T) = 0.33, 0.30, 0.27, 0.25, 0.23, 0.19 (5.1) 

for the six temperatures displayed on the bottom panel of Fig. (2). From Fig. (3), we 
see that at the last two temperatures, the Yang-Mills F"^ correlator no longer agrees with 
the strongly coupled SYM result. Therefore we conclude that for smaller than about 

^Another possibility is that the Wilson coefficient of 6 on the RHS of the OPE changes sign when r is not 
asymptotically small. This would signal a breakdown of the perturbative series for the Wilson coefficients. 
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0.25 (i.e. A smaller than about 10), one should not expect other properties of the Yang- 
Mills plasma to coincide with those of the infinite-coupling SYM plasma. In fact it is 
somewhat surprising that for 0.25 < ag < 0.30, where the function a-^{iJ,) shows a modest 
dependence on the order in perturbation theory, the thermal correlator is so much more 
similar to the strongly coupled SYM correlator than to the weakly coupled one. This may 
be related to the fact that at finite temperature, due to infrared effects, the perturbative 
expansion parameter is g rather than ctg. Thus at an energy scale where the vacuum 
polarization effects are still well approximated by the perturbative expansion, the thermal 
physics rather has a strong coupling character. 

We now discuss whether one may use the 'empirically' observed similarity to match 
the couplings of the Yang-Mills and SYM theories. By 'matching', we mean to find a way 
of comparing the two theories in such a way that they share as many properties as possible 
at a semi-quantitative level ^. Since the 't Hooft coupling is the unique parameter of the 
SYM theory, this is the only parameter we need to fix in the comparison. To what extent 
several observables can be simultaneously made similar provides a clue as to how universal 
the properties of non-Abelian plasmas are. 

The best way to match QCD with a different theory is presumably to equate a renor- 
malized quantity such as the Debye mass [21] across the two theories. However this requires 
knowing the relation between the coupling and the Debye mass on the SYM side. A tech- 
nical obstacle to this program is that the Debye mass is independent of the coupling in the 
limit of large coupling. Other observables typically lead to the same lack of sensitivity to 
A. One then needs to know 1/A corrections to the selected observable on the SYM side, 
which leads to more involved calculations and raises questions of convergence, etc. 

A different way to match the two theories is to define a running coupling based on 
a renormalized quantity, such that the weak-coupling relation holds by definition for all 
scales. One then equates the couplings of the two theories. For example, one can define 
a Yang-Mills effective coupling from the Debye mass, A(T) = 3^ in the SU(iVc) gauge 
theory, and use that value of A in the SYM theory. A priori, when the coupling is large, 
its scheme dependence is strong. For this reason, we expect that matching the observable 
itself is the superior procedure. 

Nevertheless the non-trivial agreement of the Yang-Mills and the SYM F^-correlators 
in a range of temperatures suggests that simply using the values of A = 127rQ;s = 10 ... 12, 
where the temperature-dependent values of ag are given in Eq. (5.1), is a reasonable choice 
of coupling constant to use on the SYM side. This is a moderately large coupling constant; 
for instance, the 0(1/A^/^) correction to the A = oo shear viscosity to entropy density ratio 
ri/s is about -|-50% at this coupling [64, 65, 66]. 

6. Summary 

We have found that the gluon plasma generically screens scalar and pscudoscalar fluctua- 
tions more than the vacuum does at short distance r <^ 1/T and at long distances r 3> 1/T. 

^This notion is similar to, but distinct from the (more precisely defined) matching procedure used in 
effective field theories. 
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Near however, there is a significant range of distances of order 1 /T over which the spatial 
correlations are stronger in the plasma than in the vacuum. We interpret this fact as there 
being stronger fluctuations of wavelength 0(l/r) in the plasma than in the vacuum. As 
one increases the temperature above Tc, this effect disappears soon in the scalar channel, 
but extends to about 2Tc in the pseudoscalar case. In the latter channel, the enhancement 
of these fluctuations over those of the vacuum is about a factor two. While the pseudoscalar 
and scalar channels are expected to have similar correlation functions at very short dis- 
tances and they precisely agree in the SYM theory, the two channels look rather different 
at least up to 2Tc, according to our lattice data. The scalar correlator agrees well with 
the corresponding correlator in the strongly coupled SYM theory in the range of temper- 
atures 1.2 < T/Tc < 1.9, while the pseudoscalar correlator is notably different due to the 
aforementioned strong fluctuations. These observations constitute our main results. The 
scalar fluctuations of wavelength ~ 1 /T are suppressed compared to the vacuum, while the 
pseudoscalar fluctuations are significantly enhanced. It would be interesting to see whether 
a next-to-leading order perturbative calculation would agree significantly better with the 
lattice data than the treelevel calculation does. 

We note that studying the vacuum subtracted correlators of gauge invariant operators 
at distances short compared to l/T is morally equivalent, via the operator-product expan- 
sion, to investigating the thermal expectation values of higher-dimensional renormalized 
operators. Our study thus has goals in common with the investigation of twist-two operator 
expectation values [13]. 

The semi-quantitative agreement of the scalar correlators between the pure Yang-Mills 
and the SYM theories, while the pseudoscalar channel is markedly different, highlights the 
fact that different plasmas can exhibit quite similar properties in some channels while 
differing substantially in others. 

In spite of having a reduced topological susceptibility [62] , the deconfined phase close 
to Tc exhibits strong correlations of E ■ B over distances of order l/T, which are stronger 
than in the vacuum by about a factor two. It would be interesting to see whether models 
of QCD can account for this effect. It would also be worth investigating how much this 
effect depends on the weakness of the first-order deconfining phase transition, and whether 
the effect persists at larger values of Nc, where the transition is strongly first order [67, 68]. 
A plausible mechanism for the observed strong spatial correlations is that fluctuations 
of E ■ B with a coherence length of at least l/T occur in the plasma. Perhaps these 
fluctuations have been seen in [62], where the topological lump size was found to be peaked 
at /9 ~ 1.7/Tc. This large size led the authors to conclude that this peak lies outside the 
range of applicability of their semiclassical methods. It is worth thinking about possible 
phenomenological implications of these large-amplitude, long-wavelength fluctuations, since 
the charge-separation effects of a non-zero E ■ B field configuration in the context of heavy 
ion collisions have recently received a lot of attention [69, 70, 71]. 
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A. Details of Numerical Integration 

We must numerically integrate the flow equation (3.8) 



from the horizon at r = 1 up to the AdS boundary at r ^ oo, which the initial condition 
x(r = 1) = 1. In practice we integrate to r = 20000 and verify that further increasing 
the integration domain does not change the answer. Note that while Im(x) contains a 
divergence as r ^ oo, this is a standard UV divergence that contributes only a contact 
term, and can be removed by holographic renormalization. It will not concern us, as the 
real part of x (and thus the imaginary part of Gr) has a finite limit as r — > oo. 

We cannot begin our integration at precisely r = 1 as the equation is singular there. 
We thus build a series expansion of x about r = 1: 



Plugging this into Eq. (A.l) we can determine the expansion coefficients Xn up to any 
desired order. The expressions are lengthy but straightforward to obtain and so we do not 
present them here; however we use the first three terms in this expansion to find the value 
of x{r = 1 + 6) and use this the initial condition to begin our integration at a small finite 
value of 6 {6 = 0.01 in practice). 

Note also that if we keep k finite and take a> — 0, the flow equation appears singular. 
However we know that at vanishing chemical potential the spectral density of a bosonic 
operator must be an odd function of u), and thus vanishes as a) ^ 0, although the precise 
point Co = presents numerical difficulties. Thus our code simply sets Im(Gj?(a; = 0)) = 
by hand. 
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Figure 1: Zero-temperature correlator of the trace anomaly 6{x) (top) and of the topological 
charge density q(x) (bottom). The overall normalization in the latter case is fixed by the data of 
Horvath et al. [53]. The dotted lines at small r correspond to the one- loop result with choices of 
renormalization scale (from top to bottom) /i 
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Figure 2: Thermal part of the spatial correlator of the trace anomaly 9, at Nr — 6 and sev- 
eral temperatures. The curve in the lower panel is the one-loop result (Eq. 2.7) with choice of 
renormalization scale described in section (4.3). 
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Table 2: Bare vacuum correlators for two different lattice spacings. 
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Table 3: Bare finite-temperature correlation functions at /? = 6.2822 {xs = 1.220, Zg = 1.546, 
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Table 4: Bare finite-temperature correlation functions at /? = 6.408 [xs = 1.148, Zq = 1.495, 
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